Evolution of shell structure in neutron-rich calcium isotopes 
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We employ interactions from chiral effective field theory and compute the binding energies and 
low-lying excitations of calcium isotopes with the coupled-cluster method. Effects of three-nucleon 
forces are included phenomenologically as in-medium two-nucleon interactions, and the coupling to 
the particle continuum is taken into account using a Berggren basis. The computed ground-state 
energies and the low-lying J" = 2 + states for the isotopes 42 ' 48 ' 50 ' 52 Ca are in good agreement with 
data, and we predict the excitation energy of the first J n = 2 + state in 54 Ca at f .9 MeV, displaying 
only a weak sub-shell closure. In the odd-mass nuclei 53 > 55 > 61 c a we nn d that the positive parity 
states deviate strongly from the naive shell model. 

PACS numbers: 21.10.-k, 21.30.-x, 21.60.-n, 27.40.+Z, 27.50,+c 



Introduction. - The shell model is the paradigm for 
our understanding of atomic nuclei. Doubly magic nuclei 
(i.e., nuclei that exhibit an enhanced stability) are its 
cornerstones, and they determine the properties of nuclei 
in entire regions of the nuclear chart. The magic num- 
bers - established ad hoc via a mean field plus a strong 
spin-orbit interaction more than 60 years ago by Mayer 
and Jensen for beta-stable nuclei pQ - are modified in 
neutron-rich nuclei, see for example Ref. [2] for a recent 
review. The magic nature of nuclei is reflected exper- 
imentally in enhanced neutron separation energies and 
a reduced quadrupole collectivity (i.e., a relatively high- 
lying first excited J T = 2^ state and relatively small 
electromagnetic transition probabilities from this state 
to the J x = + ground state). In doubly magic nuclei 
such as 40 Ca and 48 Ca, the J T = 2 + state appears at 
an excitation energy close to 4 MeV, while in open-shell 
calcium isotopes like 42 ' 44,46,50 Ca, this excitation energy 
is closer to 1 MeV. How these quantities evolve as we 
move towards the driplines is an open issue in ongoing 
nuclear structure research and is intimately related to our 
fundamental understanding of shell evolution in nuclei. 

For the theoretical understanding of shell evolution, 
phenomenological terms such as the tensor interaction [3] 
have been proposed. In a modern picture, three-nucleon 
forces (3NFs) play a pivotal role in shell evolution [3] . In 
the oxygen isotopes, for instance, 3NFs make 24 dou- 
bly magic and a dripline nucleus [5HZ] • Similarly, Holt et 
al. [8] showed that 3NFs are greatly responsible for the 
magic neutron number N = 28. On the other hand, ex- 
periment and theory show that the next possible magic 
number, N = 32, exhibits a smaller value for the 2 + 
excitation (but more than twice as large as seen in open- 
shell calcium isotopes) than observed in 48 Ca. This is 
often referred to as a sub-shell closure. The N = 32 
sub-shell closure is well established from experiments in 



calcium (9j [10] , titanium [IT] , and chromium [12] . How- 
ever, the situation is more complicated for neutron-rich 
calcium isotopes. For the neutron number N = 34, no 
sub-shell closure is seen experimentally in chromium |13] 
or titanium [TH [J5] , and there are some doubts regard- 
ing a sub-shell closure in calcium [16]. Different the- 
oretical predictions have been made around N = 34. 
Within the fp shell-model space, the empirical inter- 
action GXPF1 [T7j predicts a strong shell gap in 54 Ca, 
while the monopole-corrected KB3 interaction [TH] yields 
no shell gap. A low-momentum shell-model interaction 
with empirical single-particle energies and a 48 Ca core 
yields a weak sub-shell closure in 54 Ca [TS]. Shell- model 
calculations that include 3NFs predict a shell closure in 
54 Ca in the fp model space, and this shell closure is re- 
duced to a sub-shell closure (similar in strength to the 
N = 32 sub-shell closure in 52 Ca) in an enlarged model 
space that also includes the g 9 / 2 orbital [5]. Thus, the 
picture regarding the shell gap in 54 Ca is not settled yet. 
The theoretical prediction of the shell evolution in cal- 
cium isotopes is a challenging task that requires a very 
good understanding of the nuclear interaction, accurate 
treatment of many-body correlations and coupling to the 
scattering continuum j2D]- To study the shell evolution, 
we will focus on neutron separation energies, the ener- 
gies of the first excited = 2 + states, and spectra in 
the nuclei 53,55,61 Ca, which differ by one neutron from 
nuclei that exhibit a closed subshell in the naive shell 
model. 

In this Letter, we present a state-of-the-art prediction 
for the shell evolution of neutron-rich calcium isotopes. 
To this purpose, we employ nucleon-nucleon (NN) inter- 
actions from chiral effective field theory (EFT) together 
with a schematic approximation of 3NFs guided by chiral 
EFT, and utilize the coupled-cluster method to solve the 
quantum many-body problem. Chiral EFT is a system- 
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atic and model-independent approach to nuclear interac- 
tions. We employ the NN interactions at next-to-next- 
to-next-to leading order by Entem and Machleidt [2"l"]l2"2]. 
and an approximation for the chiral 3NFs that was pre- 
viously adopted in neutron-rich oxygen isotopes [7J . The 
coupled-cluster method [23] 123 is a very efficient tool for 
the computation of nuclei with a closed (sub-)shell struc- 
ture and their neighbors, and thus ideally suited for the 
task at hand. 

Hamiltonian, model space, and method. - We employ 
the intrinsic Hamiltonian 

fr_ ST ( , Mid) , M*d) \ n \ 
H - I 2mA +Vnn +y 3N 0ff ) • (1) 

l<i<j<A V 7 

Here, the intrinsic kinetic energy depends on the mass 
number A. The potential Vnn denotes the chiral 
NN interaction at next-to-next-to-next-to leading or- 
der [2"T1|2"2"] . while V3Noff is a schematic potential based on 
the in-medium chiral NN interaction by Holt et al. [21] 
(see also Ref. [21]). The potential V3Neff results from inte- 
grating one nucleoli in the leading-order chiral 3NF over 
the Fermi sphere with Fermi momentum kp in symmet- 
ric nuclear matter and is thus reminiscent of the normal- 
ordered approximation |27j . It depends formally on the 
Fermi momentum kp, the low-energy constants cp> and 
cp of the short-ranged contributions to the leading-order 
chiral 3NF, and the chiral cutoff. The latter is equal to 
the value employed in the chiral TV TV interaction [3T] . In 
the computation of neutron-rich oxygen isotopes [7] , the 
parameters kp — 1.05 fm _1 and cp = 0.71 resulted from 
adjusting the binding energies of 16 ' 22 0, while cp, = —0.2 
was kept at its value determined in light nuclei [28]. In 
this work, we use kp = 0.95 fm _1 and cp — 0.735 from 
adjusting the binding energies of 48 ' 52 Ca. It is very sat- 
isfying that the parameterization of V3Neff changes only 
little as one goes from neutron-rich isotopes of oxygen to 
the significantly heavier calcium isotopes. 

The coupled-cluster method generates a similarity- 
transformed Hamiltonian H — e~ T He T by the action 
of the cluster operator T that creates up to n-particle- 
n-hole (np-nh) excitations with respect to a reference 
state. Details of our implementation are presented in 
Refs. [2H1 ES]. We compute the ground states of the 
closed-(sub)shell nuclei 40,48,52,54,60,62 Ca in the singles 

doubles (CCSD) approximation and include n = 3 
triples perturbatively within the A-CCSD(T) approach 
of Ref. [51] . Our model space consists of up to 7V max = 19 
major spherical oscillator shells (i.e., the maximal single- 
particle excitation energy is 18 units of Tluj above the os- 
cillator ground state), and the reference state results from 
a Hartree-Fock calculation. Our basis employs oscillator 
spacings between 24 MeV < hu < 32 MeV, and in the 
largest model spaces the results we present are practically 
independent of hui. For excited states in 53 . 55 > 61 Ca above 
threshold, we use a Gamow-Hartree-Fock basis [7J [35] 



with 40 discretization points [33J. Excited states of 
the closed-shell nuclei 48 ' 52 ' 54 Ca are computed within 
the equation-of-motion (EOM) method with singles and 
doubles. The open-shell nuclei 39,4i,47,49,5i,53,55,59,6i Ca 

are computed with the particle attached/removed EOM 
methods, and we use the two-particle attached EOM 
method [3J] for the nuclei 42 < 50 < 56 Ca. Note that the em- 
ployed EOM methods are expected to reliably compute 
the separation energies and low-lying excited states as 
long as they are dominated by lp, lh, lp-lh or 2p exci- 
tations. 

Results. - Figure [I] shows the computed ground-state 
energies of the calcium isotopes and compares the results 
obtained with the Hamiltonian of Eq. |l]) to available 
data and to the results based on chiral NN interactions 
alone. The inclusion of chiral 3NFs via the in-medium 
effective potential V3Neff clearly yields a much improved 
agreement with data. The light isotopes 39 ,40,4i,42(-; a are 
slightly overbound, while the agreement is very good for 
the neutron-rich isotopes at the center of this study. The 
comparison with chiral NN forces shows that the in- 
medium effective potential V3Neff is repulsive [H] . For the 
heavier isotopes of calcium, the in-medium effective po- 
tential VaNeff becomes increasingly repulsive, and a satu- 
ration of the total binding energy sets in around 60 Ca. It 
is interesting that essentially the same interaction yields 
attraction in neutron-rich oxygen isotopes [7] . Note that 
the results for the isotopes 52_60 Ca are based on an expo- 
nential extrapolation of our results for 7V max = 14, 16, 18 
oscillator shells at the oscillator frequency Hu) = 26 MeV. 
This extrapolation adds 0.6 MeV of binding energy in 
52 Ca and 1.2 MeV in 60 Ca. For the other nuclei our 
results are practically converged with 19 shells. As a 
check, we also employed a Gamow Hartree-Fock basis for 
the s-p-d neutron partial waves and computed the ground 
states for the isotopes 54,60 Ca. The differences are small 
(up to about 0.3 MeV) and barely noticeable on the scale 
of Fig.[l] 

Figure [2] shows the energies of the J 71 " = 2 + states in 
the isotopes 42 > 48 > 50 > 52 > 54 . 56 Ca, computed in a model space 
with iV raax = 18 and Huj = 26 MeV for 48 > 52 > 54 Ca and 
N max = 16 and hw = 28 MeV for 42 ' 50 > 56 Ca. Where data 
is available, we obtain a good agreement between theory 
and experiment. The 3NFs generate the shell closure 
in 48 Ca and make N = 28 a magic number. We also 
computed the = 2 + excited state in 48 Ca using chiral 
TV TV forces alone, and we obtained E 2 + — 2.07 MeV; this 
shows that the magicity of 48 C a is due to the effects of 
3NFs in our approach. Our results predict that 54 Ca 
exhibits only a weak sub-shell closure, and this is one of 
the main results of this Letter. 

To further demonstrate this point, we also computed 
the energies of the first excited J 77 = 4 + states, and neu- 
tron separation energies S n . The results in Table [T] show 
that theory and experiment agree well for the available 
data. The theoretical results for 54 Ca suggest strongly 
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FIG. 1: (Color online) Ground-state energy of the calcium 
isotopes as a function of the mass number A. Black circles: 
experimental data; red squares: theoretical results including 
the effects of three-nucleon forces; blue diamonds: predictions 
from chiral AW forces alone. The experimental results for 
51 ' 52 Ca are from Ref. 135] . 



eral levels in 52 > 53 > 54 > 55 > 56 Q a ^ rjv^g sca ttering continuum 
is shown as a gray band. Note that we computed the 
jTc _ 2+ excited state of 54 Ca as a lp-lh excitation of 
the 54 Ca ground state and as 2p excitations of the 52 Ca 
ground state. In 17 major oscillator shells we obtain 
E 2 + =1.80 MeV and E 2 + =1.90 MeV, respectively. The 
two methods are in good agreement with each other, and 
this shows the quality of the employed methods. 
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FIG. 3: (Color online) Theoretical excitation spectra of 
52_56 Ca compared to available data. The continuum is in- 
dicated as a gray area. 



In the naive shell-model picture, the order of single- 
particle orbitals above the Fermi surface of 60 Ca is 59/2, 
ds/2 and s\/2- Mean- field calculations predict that the 
corresponding canonical orbitals are close in energy [38j . 
yielding deformation in the isotopes 60_70 Ca 39 . The 
isotope 58 Ca is the heaviest isotope of calcium that has 
been produced [40 , and 52 Ca is the most neutron-rich 
isotope whose mass is known [35]. Nuclear energy func- 
tionals and mass models [35] [35] predict the neutron 



FIG. 2: (Color online) Excitation energies of J w = 2 + states 
in the isotopes 42 > 48 > 5 °. 52 . 54 . 56 c a (experiment: black, theory: 
red). 



that this nucleus will only exhibit a weak sub-shell closure 
but not a shell closure. The neutron separation energy 
in 52 Ca is a noteworthy example. The atomic mass table 
evaluation [3S] gives the value S„ = 4.7 MeV (estimated 
from systematics) for this nucleus. The very recent mea- 
surement [35], however, yields S n « 6 MeV, in much 
better agreement with our prediction of S n = 6.59 MeV. 

To further study the evolution of shell structure we 
focus on spectra. Figure [3] shows the computed spec- 
tra for 52_56 Ca, and compares them to available data. 
In 52 ' 53 Ca, our calculations suggest spin assignments for 
measured levels [37], while we give predictions for sev- 





48 Ca 


52 Ca 


54 Ca 


E 2+ (CC) 
E 2+ (Exp) 


3.58 
3.83 


2.19 
2.56 


1.89 
n.a. 


E 4+ /E 2+ (CC) 
£ 4+ /£ 2+ (Exp) 


1.17 
1.17 


1.80 
n.a. 


2.36 
n.a. 


S n (CC) 
S*„(Exp) 


9.45 
9.95 


6.59 
6.0* 


4.59 
4.0 f 



TABLE I: Excitation energies E 2 + (in MeV) of the lowest- 
lying J n — 2 + states, neutron-separation energies S n (in 
MeV), and ratio of energies E A + / E 2 + from coupled-cluster 
theory (CC) compared to available data (Exp) [n.a. = not 
available; * = from Ref. [35] : ^=from atomic mass table evalu- 
ation [36]]. The theoretical results point to a sub-shell closure 
in 54 Ca. 
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drip line to be around 70 Ca. Figure [3j shows that the 
J 17 = 5/2 + resonant excited states in 53,55 Ca are lower 
in energy than the J 71 " = 9/2+ states, and this is a de- 
viation from the naive shell-model picture. There is no 
J n = 1 /2 + resonance in the continuum due to absence of 
a centrifugal barrier. Due to the large I = 4 centrifugal 
barrier the J n = 9/2+ resonances can be considered as 
quasi-bound states, while the J 71 ' = 5/2 + has a signifi- 
cantly larger width, see Table [TI| for details. Our results 
shown in Fig.[T]yield 60 Ca unbound with respect to 56 Ca. 
However, we cannot rule out the existence of 60 Ca since 
correlations beyond triples in the coupled-cluster expan- 
sion may play a larger role for 60 Ca than for the neighbors 
of 54 Ca. Our computations (including the scattering con- 
tinuum) suggest that 61 Ca is a very interesting nucleus. 
We find resonant states J* = 5/2+ and J 7 * = 9/2+ at 
excitation energies of 1.1 MeV and 2.2 MeV above thresh- 
old, respectively, and a virtual J 7 ' = 1/2+ state practi- 
cally at threshold. This ordering of the resonances is con- 
sistent with the results for 53,55 Ca. Table ITT1 summarizes 
our results for the J 71 " = 5/2 + and J v = 9/2+ resonances 
in 53 - 55 ' 61 Ca. We also find the 62 Ca ground state to be 
very close to threshold, unbound by about 0.2 MeV with 
respect to 60 Ca and entirely dominated by (si/ 2 ) 2 con- 
figurations. In our calculations, the correct treatment 
of the continuum is essential: (i) in 55 Ca the scattering 
continuum lowers the energy of the excited J 77 = 5/2 + 
by about 2 MeV when compared to a harmonic oscilla- 
tor basis, (ii) in 61 Ca, the oscillator basis yields the level 
ordering of the naive shell model. 





53 Ca 


55 Ca 


61 Ca 


r 


Re[E] 
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Re[E] 
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Re[E] 
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5/2+ 


1.99 


1.97 


1.63 


1.33 


1.14 


0.62 


9/2+ 


4.75 


0.28 


4.43 


0.23 


2.19 


0.02 



TABLE II: Energies of the 5/2 and 9/2 resonances in 
53,55,6iQ a R e [J5] j s the energy relative to the one-neutron 
emission threshold, and the width is T = — 2Im[_E] (in MeV). 

As a check, we also computed the lowest excited states 
with spin and parity r = 2+, 4+, 6+ in 50,54,56 Ti Thege 
states result from attaching two protons to the closed- 
core reference nuclei 48,52 ' 54 Ca, respectively, and we em- 
ploy 2p-Qh and Zp-lh excitations in the computation. 
Figure [4] shows a reasonably good agreement with the 
data, with a maximal deviation of about 0.6 MeV for the 
2\ state in 56 Ti. Let us compare the titanium isotopes to 
the isotopes 50 > 54 > 56 Ca which might also be computed by 
attaching two neutrons to 48 ' 52 ' 54 Ca, respectively. The 
analysis of the cluster amplitudes shows that the 2p-0h 
cluster amplitudes account for about 60% of the total am- 
plitude for the isotopes of titanium, while this number is 
about 72% for the calcium isotopes (with the remain- 
ing weight carried by 3p-lh amplitudes). It is insight- 
ful to compare the ratios of 3p-2h to 2p-0h amplitudes. 



These are about 0.67 in titaniums and only 0.39 in calci- 
ums. Thus, the isotopes of titanium are more correlated 
due to proton-neutron interactions than the correspond- 
ing calcium isotopes. We believe that one needs 4p-2h 
amplitudes for a more accurate description. 
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FIG. 4: (Color online) Theoretical excitation spectra of 
50,54,56 Ti compared to data 

Summary. - We employed interactions from chiral ef- 
fective field theory and performed coupled-cluster com- 
putations of neutron-rich calcium isotopes. The evolu- 
tion of shell structure is understood by our computations 
of neutron-separation energies, low- lying J 77 — 2+,4+ 
states, and excitations in the odd-mass nuclei 53,55 ' 61 Ca. 
Our results confirm that the shell closure in 48 Ca is due 
to the effects of three-nucleon forces, and that 52 Ca ex- 
hibits a sub-shell closure while 54 Ca exhibits only a weak 
sub-shell closure. We make several predictions for spin 
and parity assignments in the spectra of very neutron- 
rich isotopes of calcium and predict a level ordering in 
53,55,61 Q a -|-h a t j g a £ var j ance w ith the naive shell model. 
Overall, we find that effects of three-nucleon forces and 
the scattering continuum are essential in understanding 
the evolution of shell structure towards the dripline. 
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